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Abstract 

We show that the rather unexpected pressure dependence of superfluid density observed near 
the superfluid-supersolid transition by Kim e£.a/.[M.H.W. Chan, private communication], can be 
understood if the transition from superfluid to supersolid state is a second order or weakly first order 
transition from the superfluid state to a super-CDW state with non-uniform Bose-condensation 
amplitude. The suppression of superfluid density is a direct consequence of softening of phonon 
mode at finite wave- vector \Q\ ~ Qq around the quantum phase transition. 



Following the basic theoretical understanding of superfluidity in in ife 4 , the possibility 
of a so-called supersolid state has been propose d 1 ! 2 ' 3 ' 4 . In its simplest form, the theoretically 
proposed supersolid state can be thought of as a state with spatially non-uniform Bose- 
condensation amplitude ^(af) 1 ! 2 ' 3 ' 4 . For example, 



(60 > bg), corresponding to a Bose-condensed charge-density wave, or super-CDW, state. 
Notice that this state, if stable, has a non-zero superfluid density by construction, and is 
very different from the usual insulator state where superfluid density vanishes. 

The supersolid state is observed recently in a torsional oscillator experiment by Kim 
et.a l 11 ' * 12 where finite (but small) superfluid density was observed in an otherwise solid He 4 
state. A rather unexpected finding of the experiment is that the superfluid density is not 
a monotonic function of pressure but increases first as pressure increases, and decreases 
again when pressure is higher, whereas naively we expect that the superfluid density should 
be a monotonic decreasing function with pressure 13 . A number of theoretical works were 
triggered by the discovery of this new stat o 5 i 6 i 7 i 8 i 9 i 10 . The purpose of our paper is to point 
out that the rather unexpected pressure dependence of superfluid density as found in the 
experiment can be understood if we identify the transition from superfluid to supersolid 
state as a second order or weakly first order phase transition from the superfluid state to 
a super-CDW state, whereas the transition to the real insulator state occurs at a higher 
pressure. The suppression of the superfluid density p s is a direct consequence of enhanced 
thermal fluctuation coming from softening of phonon mode at finite wave- vector \Q\ ~ Q 
around the quantum phase transition and is quite independent of the microscopic details of 
the system. 

To demonstrate we consider an effective action for the boson system of form 



where p = p + p\ and 9 are the density and phase variables characterizing the boson system 
and V[p] is an energy functional. The ground state density profile p of the system is deter- 
mined by minimizing V[p] with respect to p. pi describes the density fluctuations above the 
ground state. The boson creation (annihilation) operators are related to the density-phase 



if)(x) = b + b Q cos(Q.x) 



(1) 
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variables by ip + {ip){x) = J p{x)e W^W. Specifically, 

V[ P }= Jd 3 x^-(VVP) 2 + U[pl (3) 

in usual interacting boson problem where U[p] — hj d 3 x J d 3 x'p(x)U(\x — x'\)p{x') is the 
potential energy of the bosons with U(x) being the corresponding interaction potential. 
Here we shall treat V[p] as a phenomenological effective potential with which the ground 
state and low energy properties of the system can be described by the Gaussian theory. 

In the normal superffuid state, the density is uniform and p(x) = n. Small density 
fluctuations around the ground state can be described by expanding V[p = n + Ap] in a 
power series of Ap^, i.e. 

V[p] ~ V[n]+Y,A(\Q\)Ap Ap_$+ £ Q^A^Ap^Ap^^ (4) 

Q Qi,Q2 

+ E ^ C{Q u Q 2 ,Qz)A PQ Ap3Ap Q Ap_ 0i _$ 2 _Q 3 + .... 

Qi,Q2,Q3 

where we assume that the system is translational invariant. Notice that the coefficients 
A,B,C,... are in general functions of external parameters like pressure and temperature. 

For B = 0, the transition from the normal superffuid state to super-CDW state is charac- 
terized by a change of sign of A{Q) from positive to negative at some wavevector Q = Q > 0. 
The change in A(Q) can be induced by pressure as in the experiment of Chan eta/ . 11 ' 12 or 
through other means. In this case the quantum phase transition is a second order phase 
transition. For B ^ 0, the transition becomes a first order phase transition occurring at 
finite value of A(Q ) > 0. In the following, we shall show that complete suppression of 
superffuid density at the quantum phase transition is a natural consequence of the above 
phenomenological Lagrangian when B = 0, and suppression of superffuid density still occurs 
when B is small, when the transition is weakly first order. Notice that for the super-CDW 
state to be stable, A(Q) must has the property that it is less than zero only at a small region 
around Q ~ Q , and is positive otherwise. 

To get some physical feeling of A(Q), we consider V[p] of form (JHJ). In this case, it is 
easy to show that A(Q) = (e^ + 2nU(Q))/An, where = h 2 q 2 /2m is the kinetic energy 
and U(q) = J d d xe^ x U(x) is the Fourier transform of the inter-particle interaction U(x). 
In terms of the inter-particle potential U(q), the stability requirement of the CDW state 
means that we must have U(q —>■ 0) > with an absolute minimum at q = Qq ~ (particle 
spacing) -1 , with e^ Q + 2nU(Q ) < around this region. 
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It is interesting to note that this form of effective interaction potential has been used 
to describe low temperature properties of helium liquid^ and can be derived from, for 
example, the STLS scheme^. In these theories, U(Qq) < leads to a peak in the structure 
factor S(q ~ Qo) and a corresponding dip in the phonon spectrum E(q ~ Qq)- As pressure 
increases, U(Q ) decreases further and the peak magnitude in S(Q ) increases. The increase 
in magnitude of S(Qo) is interpreted as a precursor to formation of Wigner crystal^^. 
Notice that correspondingly the phonon energy E(q ~ Qo) softens. The superfluid-super- 
CDW instability is driven by an instability in the phonon spectrum at q ~ Qq when E(Q ) — > 
and S(Qo) diverges. 

We next analyze the phenomenological Lagrangian (j2J) and (@J) in Gaussian approxima- 
tion. First we consider B = and A(Q) > 0. In this case higher order terms in (@J) are 
unimportant. After integrating out the pi (or Ap) fields, we obtain 



f P dr J d d q LeMo)) 2 + wr^m) 2 ) • (5) 



So 



where the usual phonon dispersion E(q) = Jne^A(q) ~ J €g> (e^ + 2nU (q)) follows. As 



q — > 0, E(q) — > ynA(0)e^ ~ which is the usual superfluid Goldstone mode. We also 
observe that E(Qq — > 0) — > when A(Qq) —>■ 0, corresponding to softening of the phonon 
spectrum at the \Q\ = Qo when the CDW instability occurs. 

It is straightforward to compute the depletion in superfluid density near the superfluid- 
super-CDW transition coming from thermal fluctuations. The superfluid density is given by 
p s (T) = n - p n (T), where 

Pn(T) = "o / 7^(^ —J— U=E(0, (6) 



3 J (2tt) 3V ^ de 

E 

is the normal liquid density 17 . ub{E) = [e k B T — is the boson occupation number. 

To proceed further we assume E(q) 2 ~ E(Q ) 2 + c 2 (q — Qq) 2 around q ~ Q , where 
c 2 = d 2 E(q)/dq 2 \ q= Q > and E(Q )(A(Q )) — > 0. The contribution to p n from around 
\Q\ ~ Qo can be estimated by evaluating the integral. For k B T >> E(Q ) we obtain 

where p^ h (T) ~ T 4 is the contribution from ordinary phonon around zero momentum and 



^^m^^, (8 ) 

67rc V 7rk B T 



for ksT « E(Qo). We observe that p n (T) diverges as IcbT/ E(Q Q ) at the superfluid-super- 
CDW transition when ksT » E(Q ) — > 0, indicating that superfluidity is destroyed as a 
result of thermal fluctuations at the critical point at any non-zero temperature T ^ 0. The 
suppression of superfluidity is analogous to the suppression of Bose-condensation by thermal 
fluctuation in one- dimensional Bose gas. Notice that the suppression occurs because the 
instability occurs at finite wave-vector Qo and would be absence in transitions where the 
instability occurs at Q = or in lattice models where Q commensurate with the underlying 
lattice structure. 

The ksT / E(Qq) singularity is cutoff by a first-order phase transition when B ^ 0. To 
see how this occurs we consider the potential energy (@J), assuming that A(Q) approaches 
zero at Q = Q and the boson system forms a super-CDW state with period ~ Qq 1 - At 
the vicinity of the transition the CDW wave can be viewed as superposition of density 
waves with same amplitude but travelling at different directions with different phases^, i.e., 
p(x) ~ So Ap^e lG ' x where G's are primary reciprocal lattice vectors with |G| = Qo and 
Apg ~ \Ap \e w G. In this case, the potential energy (jlj becomes 



state, B ~ J2g's B (^i, G 2 )e i ^ +e ^ and C ~ £^, s C (^i, G 2 , G 3 )e ii9l+d2+9 ^ > are constants 



that depends on Q , M and the lattice structure of the CDW state. 

Minimizing V[p] up to fourth order terms, it is easy to see that Ap ~ B/2C > when 
MA < and the second order phase transition is cutoff by a first order phase transition 
when B ^ 0. Correspondingly, the fc^T/ E(Q ) divergence of p n (T) computed in Eq. is 
cutoff with 



Notice that the superfluid density changes discontinuously across the phase transition as 
a result of the discontinuous change in the excitation spectrum when B ^ 0. The superfluid 
density in the super-CDW state can be written as p GDW (T) = p GDW (0) — p GDW (T), where 
p c DW (0) ~ n - a(Apo) 2 and p% DW {T) ~ p n {T) + (3{T){Ap ) 2 around the phase transition, 
where a,(3(T) are non-diverging constants depending on the CDW lattice structure. It is 
easy to show that at the CDW side of the transition, V[p] ~ V[n] + V[Ap ] + A'pf, where 



V[p] - V[n\ + MA(Q )(Ap ) 2 + B(Q )(Ap ) 3 + C(Q )(Ap ) 4 + .... 



(9) 



— * 

where M = number of reciprocal lattice vectors G's that enters the construction of the CDW 
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A' = B 2 /AC. Therefore we expect that the low energy cutoff is replaced by 

Ecdw(Qo) 

in this case. Phonon excitations with energy above E C dw(Qo) is on ly weakly affected by the 
presence of Ap ^ and p n (T) is given by Eq. ((7J) and except that E(Q ) is replaced 
by E CDW (Q Q ). 

When applying our results to the superfluid-supersolid transition, it suggests that what 
Kim et. al-~& have observed experimentally may be the depletion of superfluid density around 
the quantum phase transition due to thermal fluctuations. Moreover, it suggests that similar 
depletion will be found in the superfluid side of the transition. The depletion is temperature 
dependent and there will be a discontinuity in the superfluid density across the transition 
if the transition is first order, with the superfluid density in the normal superfluid side be- 
ing suppressed more, at least when B is small enough. The phonon dispersion at Q ~ Qo 
hardens*'- and the superfluid density rises again when the system moves away from both sides 
of the critical point. Notice equations (7) and (8) predict also a rather specific tempera- 
ture dependence of superfluid density around the transition. Therefore, a measurement of 
temperature-dependent superfluid density across the transition can determine the order of 
the transition, and provides a test to our theory. 
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